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Abstract. A complete classification of quaternary complex Hadamard matrices of orders 
10, 12 and 14 is given, and a new parametrization scheme for obtaining new examples of 
affine parametric families of complex Hadamard matrices is provided. On the one hand, it 
is proven that all 10 x 10 and 12 x 12 quaternary complex Hadamard matrices belong to 
some parametric family, but on the other hand, it is shown by exhibiting an isolated 14 x 14 
matrix that there cannot be a general method for introducing parameters into these types 
of matrices. 
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1. Introduction 

A complex Hadamard matrix H of order n is an n x n matrix with unimodular entries 
satisfying HH* = nl where * is the conjugate transpose, / is the identity matrix of order 
n, and where unimodular means that the entries are complex numbers that lie on the unit 
circle. In other words, any two distinct rows (or columns) of H are complex orthogonal. 
Complex Hadamard matrices have various applications in mathematics ranging from cod- 
ing and operator theory [12] to harmonic analysis P]. They also play a crucial role in 
quantum information theory, for construction of quantum teleportation and dense coding 
schemes [19], and they are strongly related to mutually unbiased bases (MUBs) [21]. In this 
paper we are primarily concerned with n x n Butson-type Hadamard matrices [2], denoted 
by BH(g,n), which are complex Hadamard matrices composed of gth roots of unity. The 
notations BH(2,n) and BH(4,n) correspond to the real Hadamard matrices and quaternary 
complex Hadamard matrices, respectively. 

Recently there has been a renewed interest in BH(g, n) matrices. In particular, they 
have been used as starting-point matrices for constructing parametric families of complex 
Hadamard matrices for various q and n. In a series of papers Di^a [4], SzoUosi [E], and 
Tadej and Zyczkowski [V7\ obtained new, previously unknown parametric families of complex 
Hadamard matrices from BH(g, n) matrices. The parametrization of a BH(g, n) matrix is 
an operation where some of the matrix entries are replaced by infinitely differentiable, or 
smooth, functions mapping a real-valued or complex-valued argument vector into the set of 
unimodular complex numbers. Any assignment of values to the parameters yields a complex 
Hadamard matrix and some assignments produce BH(g, n) matrices. The parametrization of 
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the BH(g, n) matrices makes it possible to escape equivalence classes and therefore to collect 
many inequivalent matrices into a single parametric family. Also the switching operation 
[16] . a well-known technique in design theory, has this property. 

The existence of BH(g, n) matrices is wide open in general. Even the simplest case for 
g = 2 is undecided, as the famous Hadamard conjecture, stating that BH(2, 4k) matrices exist 
for every positive integer k, remains elusive despite continuous efforts [3J. Real Hadamard 
matrices, or BH(2,n) matrices, are completely classified up to = 28, and there has been 
some promising advance in enumerating the case n = 32 very recently (see [7] and the 
references therein). For other values of q we have some constructions [2] and some non- 
existence results [20] ■ Harada, Lam and Tonchev classified all 16 x 16 generalized Hadamard 
matrices over groups of order 4 and obtained new examples of BH(4, 16) matrices This 
particular result motivated the authors to investigate the existence of BH(4, n) matrices in 
more general and to start enumerating and classifying them for small n. The census of 
quaternary complex Hadamard matrices up to order 8 were carried out in [H]. The aim 
of this paper is to continue that work, and give a complete classification of the BH(4, n) 
matrices up to orders n = 14. 

This work has two major parts: first, we classify all BH(g, n) matrices using computer- 
aided methods, and secondly, we define parametric families of complex Hadamard matrices 
by introducing parameters to the BH(g, n) matrices. Parametric families offer a compact 
way of representing a large number of BH(g, n) matrices, and they also yield information 
about complex Hadamard matrices. 

In the previous work [13] the authors collected all BH(4, n) matrices contained in various 
parametric families from the existing literature and then confirmed with an exhaustive com- 
puter search that those are the only examples. Here the situation is exactly the opposite, as 
it turns out that almost all BH(4, tt.) matrices of orders n = 12 and 14 found in the current 
work by the computer search are previously unknown. Therefore we are facing the inverse 
problem: we need to encode a given collection of BH(4, n) matrices by parametric families. 

The outline of the paper is as follows. In Section |2] we give some basic definitions and 
results that are used in later sections. Then in Section [3] we recall the computer-aided 
classification method of difference matrices from [10] and highlight the main differences and 
necessary modifications required for our purposes. In Section H] we introduce a new method 
for parametrizing complex Hadamard matrices and recall some relevant materials from the 
existing literature. Finally, in Sections [5l |6l and [7] we present the classification of BH(4, n) 
matrices for n = 10, 12, and 14, respectively. To improve the readability of our paper a rather 
long list of complex Hadamard matrices has been moved from the main text to Appendix |X1 
The matrices the authors obtained are also available in an electronic format in the on-line 
repository [9] as a supplement. Throughout this paper we adopt the notations from [T7] for 
known complex Hadamard matrices, such as D\q {a,b,c), etc. 

2. Preliminaries 

Two BH(g, n) matrices are equivalent if the first matrix can be transformed into the second 
one by permuting the order of its rows and columns and multiplying a row or a column by 
some gth root of unity. The automorphism group of a BH(g, n) matrix H is the group of 
pairs of monomial matrices (P, Q) such that H = PHQ; a monomial matrix is an n x n 
matrix having a single nonzero entry in each row and column, these nonzero entries being 
complex qth roots of unity. Note that the automorphism group depends on the choice of q. 
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We employ invariants for finding the equivalence class of a BH(g, n) matrix in a computa- 
tionally efficient way. The determinant of a kxk submatrix of matrix H is called a vanishing 
minor, or a zero minor, of H if the determinant is zero. Throughout this paper we shall 
repeatedly use the following 

Lemma 2.1. Let H be a complex Hadamard matrix of order n and 2<k<n — 2be 
an integral number. Then the number of k x k vanishing minors of H is invariant, up to 
equivalence. 

The number oik x k vanishing minors is just a special case of a more powerful invariant, 
the fingerprint [15], but it is sufficient for our purposes. 

The g-rank, or Z^-rank, of an integral matrix L of order n is the smallest positive integer r 
such that there are integral matrices S and T of orders nxr and rxn respectively, such that 
ST = L (mod q). The Z^-rank of a BH(g, n) matrix H is the Z^-rank of the (0, 1, . . . , g — 1)- 



matrix L for which H = EXP ( —L ) , where EXP denotes the entry-wise exponentiation 



function. It is easy to see that the Zg-rank is, again, invariant up to equivalence. A class of 
BH(g, n) matrices with small Zg-rank has some interesting applications in harmonic analysis 



As it seems that a complex Hadamard matrix H shares every important property with its 
hermitian H*, complex conjugate H and transpose H^, we introduce the concept of ACT- 
equivalence. Two complex Hadamard matrices H and K are called ACT-equivalent, if H 
is equivalent to at least one of K, K*, K or . The concept of this refined equivalence 
simplifies the presentation of our results as we can avoid unnecessary repetitions in our 
summarizing tables [9]. 



The computer-aided methods employed in the classification of Butson-type Hadamard 
matrices in this work are very similar to the methods used for the classification of difference 
matrices over cyclic groups [10]. Therefore, we give here only a summary of the relevant 
ideas and describe in detail only the points where this work differs from the work done on 
difference matrices. 

We perform an exhaustive computer search of the space of all Butson-type Hadamard 
matrices BH(g, n) with the given parameters g = 4 and n = 10, 12, 14. Because the size of 
the search space grows very quickly as the parameters are increased, we have to prune most 
of the search space and we do this by considering only inequivalent matrices. 

An mxn matrix over gth roots of unity is said to be a candidate BH(g, n) matrix, denoted 
by BH(g, m, n), if the rows are orthogonal as in a BH(g, n) matrix, but the number of rows 
m < n. A BH(g, m, n) matrix is dephased (or normalized) if the first row and the first 
column contain only the value 1. Since every BH(g, m, n) matrix is equivalent to a dephased 
one, only these matrices are considered in the computer search. 

The dephased candidate BH(g, m, n) matrices are organized into a tree where the child 
nodes of a BH(g, m, n) are all BH(g, m + 1, n) that are obtained by appending a row to the 
node, cf. [TUl Definition 2.6]. An exhaustive search of this tree is performed by weak canonical 
augmentation [6], Section 4.2.3], which has the advantage over the simpler breadth-first search 
used in [lOj that it can be performed easily in parallel. 




Consider a tree of dephased candidate BH(g, m,?T,) matrices. Let p{X) denote the parent 
of the node X in the tree. Every node X in the tree has a finite sequence of ancestors from 




3. Constructing Butson-type Hadamard matrices 
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which it has been constructed by appending a row to it: 

(1) X,piX),p{piX)),pipipiX))),... 

In the tree each X occurs only once but there are typically many nodes Y with X = Y. 
Such a node Y has also a sequence of ancestors from which it has been constructed: 

(2) Y,piY),pipiY)),pipipiY))),... 

Even though X = Y the ancestor sequences ([T]) and ([2]) need not consist of the same nodes 
up to equivalence. This means that the sequences ([T]) and ([2]) can be distinct on the level of 
equivalence classes even if X = Y. The main idea is now to exploit such differences among 
equivalent nodes in rejecting equivalent nodes. 

Let Tnr denote the set of all non-root nodes in the search tree T. Associate with every 
object X G Tnr a weak canonical parent w{X) G T such that the following property holds: 

(3) for all X, Ye Tnr it holds that X = Y implies w{X) = w{Y). 

The function w defines for every non-root node X a sequence of objects analogous to ([T]): 

(4) X, w{X),w{w{X)),w{w{w{X))), . . . 

Because of ([3]), any two equivalent objects have identical sequences (jl]) on the level of 
equivalence classes of matrices. When the search tree is traversed in depth-first order, a 
node X and the subtree rooted at it is considered only if it has been constructed in the 
canonical way specified by (jl]); that is, every matrix in the ancestor sequence ([1]) should be 
equivalent to the matrix in the corresponding position in the sequence (jl]). By (j3]) we obtain 

p{X) = w{X) ^w{Y)=p{Y) 

In other words, equivalent matrices generated by weak canonical augmentation have equiv- 
alent parent matrices, and assuming that the same holds for the parents, this implies that 
equivalent matrices must be siblings in the search tree. This reduces the size of the search 
tree dramatically. 

The problem of checking equivalence of BH(g, m,r;,) matrices is solved by transforming 
it into a corresponding graph isomorphism problem in exactly the same way as was done 
with difference matrices over cyclic groups in [lOJ, Section 3]. Each BH(g, m,n) matrix is 
mapped to a directed graph, called the equivalence graph of the matrix. Two BH(g, m,n) 
matrices are equivalent if and only if their equivalence graphs are isomorphic graphs. The 
definition of an equivalence graph is obtained from the definition of a difference matrix graph 
[To] Definition 3.2] by considering gth roots of unity as a cyclic group where the operation 
is complex multiplication. 

The software package Nauty [11] is used for checking the isomorphism of equivalence 
graphs. For each equivalence graph Nauty calculates a graph, called the canonical graph, 
which has the property that two equivalence graphs are isomorphic if and only if they have 
the same canonical graph. This implies that two BII(g, m, n) matrices are equivalent if and 
only if the canonical graphs of the equivalence graphs of the matrices are the same. We 
define the canonical graph of a matrix X as the canonical graph of the equivalence graph of 
X. 

Transforming a matrix to a graph yields also a weak canonical function w having the 
property (j3]). Let P{X) be the set of matrices obtained by removing a row from a matrix 
X G Tnr and let <g denote a total order on the set of canonical graphs. We define w{X) as 
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the matrix Y G P{X) which has the smallest canonical graph under <g among the canonical 
graphs of matrices in P{X). 

The definition of equivalence of BH(g, n) matrices given in Section |2] agrees with the 
equivalence relation =* defined in [TO', Section 2] and the equivalence inducing group E* in 
[ini Section 5]. 

To get confidence in the computational results, we perform a consistency check by counting 
the number of candidate and complete BH(g, n) matrices in two different ways. In the 
first method the sizes of the equivalence classes of BH(g, m,?T,) matrices are obtained from 
inequivalent BH(g, m, n) matrices and the automorphism groups of their canonical graphs 
via the orbit-stabilizer theorem. The second method determines the number of BH(g, m, n) 
matrices from the number of BH(g, m — l,n) matrices and number of ways these smaller 
matrices can be augmented to yield a BH(g, m,n) matrix. This double-counting method 
is the same as in [lOl Section 5] except that the analytical formula for 2 x c difference 
matrices and the reduction used for 3 x c difference matrices are not valid for BH(g,2,n) 
and BH((3', 3, n) matrices, respectively. With BH(g,n) matrices the double-counting starts 
one level earlier with BH((3', l,n) matrices, which all belong to the same equivalent class of 
size g". 

The exhaustive computer search yields all inequivalent BH(g, n) matrices and the auto- 
morphism group for each matrix. 

4. Parametrizing complex Hadamard matrices 

Once a complete set of inequivalent BH(4, n) matrices were obtained the authors inves- 
tigated whether these matrices lead to parametric families of complex Hadamard matrices. 
We recall the following basic method to introduce free parameters into complex Hadamard 
matrices of even order |15l Lemma 3.4]. 

Lemma 4.1. Let H he a dephased eomplex Hadamard matrix of order n and suppose that 
there exist a pair of rows in H , say u and v, such that for every i = 1,2, . . . ,n, uf = vf . 
Then for all such i for which Ui + Vi = replace Ui with aui and Vi with avi, where a is a 
unimodular complex number to obtain a one-parameter family of complex Hadamard matrices 
H{a). 

Lemma 14.11 is a general method for introducing affine parameters into real Hadamard 
matrices, irrespectively of their order, and it can be applied to various other matrices as 
well. Unfortunately, however, it cannot be applied to BH(4, n) matrices when n = 2 (mod 
4). In the following we present a method working for some of these BH(4,n) matrices. 

Let us denote by 1"* the all-1 row vector of length m, and let (., .) denote the standard 
inner product in with the convention that it is linear in the first, and conjugate linear in 
its second argument. 

Theorem 4.2. Let H be a dephased complex Hadamard matrix with the following block 
structure 



1 


1 


1 


IP 




1 


a 


b 


X 


y 


1 


b 


a 


X 


-y 


vY 




T 

Z 


A 


B 


T 

W 


T 

—w 


C 


D 



6 PEKKA H.J. LAMPIO, FERENC SZOLLOSI AND PATRIC R.J. OSTERGARD 

where a and b are arbitrary unimodular numbers. Then, after replacing the row vectors y 
with ay and w with aw we obtain a one-parameter family of complex Hadamard matrices 
H{a) where a is unimodular. If, in addition, b = a we can continue by replacing w in H{a) 
with a(5w to obtain a two-parameter family of complex Hadamard matrices H{a,P) where a 
and (5 are unimodular. 

Proof. We need to show that the rows of H{a., /3) are pairwise orthogonal. From the orthog- 
onahty of the first three rows of H we get 

([l,l,l,lM''],[l,a,6,a;,y]) = 1 

([l,l,l,F,l'?],[l,6,a,a;,-y]) = oj ^ ^''^^ ^' 

and 

= ([1, a, 6, X, y], [1, 6, a, X, — y]) = 1 + ah + ba + p — q 

= ([l,a,6,x,a?/], [l,a,6,x, -ay]) , 

and hence the first three rows of H[a,f3) are pairwise orthogonal. Similarly, it is easily seen 
that the rest of the rows (beyond the first three) are pairwise orthogonal within themselves. 
Additionally, the first row is trivially orthogonal to all further rows. Therefore it remains to 
be seen that the second and third rows of H{a,P) are orthogonal to all rows below them. 

We show first that they are orthogonal to the rows which are of type [1, Zi, Zi, Ai, Bi], 
i = 1, . . . ,r. In the original matrix H (i.e., prior to parametrizing) we have 

(5) 1 + Zi{a + b) + {Ai,x) + {Bi,y) = 0, 



(6) l + z,(a + b) + {Ai,x)-{Bi,y)=0, 

and hence {Bi, y) = for every i = 1, . . . ,r. It follows that after parametrization equations 
(|5]) and ([6]) remain valid. 

We proceed by proving that rows that are of type [1, Wi, —Wi, Ci, Di], i = 1, . . . ,s, after 
parametrization, are orthogonal to the second and third row of H(a,(3). Again, in the 
original matrix H we have 

(7) 1 + Wia-Wib+ {Ci,x) + {Di,y) =0, 



(8) l-Wia + Wib+{Ci,x)-{Di,y) = 0, 

and hence {Ci, x) = —1 for every i = 1, . . . , s. It follows, that ([7]) are (|8]) are valid, provided 
that 

(9) Wia-Wib+{Di,y) = 

holds and therefore ([7]) and ([8]) remains true, after parametrization. If, in addition, b = a, 
then {Di,y) = for every i = l,...,s, and hence ([9]) holds, independently of the scalar 
factor in w. 

The one-parameter family H{a) can be considered as H{a,a). The equations (E]) - ([H]) 
hold as the condition a = 6 is not required for them. From the original matrix H we get 

WiO — Wib + {Di, y) = =^ aWiO — awib + {Di, ay) = 0, 

and holds for H{a,a) also when a ^ b. □ 
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Remark 4.3. If a and b are as in Theorem 14.21 then it is easy to see that the real part of ab is 
an integral number, and therefore b G ±a -{1,0;, w^, i}, where a; is a primitive complex third 
root of unity. Therefore one hopes to apply the parametrizing scheme described for complex 
Hadamard matrices with fourth and/or sixth roots of unity. 

Theorem 14. 21 describes a local property of the complex Hadamard matrix H. Its conditions 
can be fairly easily checked, even by hand, and it can be implemented as a computer program 
to construct infinite families automatically. 

It is natural to ask how many degrees of freedom (i.e., independent, smooth parameters) 
can be introduced into a given complex Hadamard matrix. We recall the following funda- 
mental concept from [T7j: The defect d{H) of an n x n complex Hadamard matrix H reads 
d{H) = m — 2n + 1, where m is the dimension of the solution space of the following real 
linear system with respect to the matrix variable R G R"'^'^: 

n 

(10) H^,kHj,k {Ri,k - Rj,k) = 0, l<i<j<n. 

k=l 

In what follows we motivate the formula (fTOj) . Consider a dephased complex Hadamard 
matrix H and a phasing matrix R whose first row and column is 0, and all other entries are 
a real variables. Now consider the matrix K := H o EXP (i-R), where o is the entry- wise 
product and EXP is the entry-wise exponential function, respectively. The resulting matrix 
K is the most general parametrized matrix, stemming from H. Note that K is unimodular, 
but not necessarily complex Hadamard. In order to ensure this latter condition we spell out 
the orthogonality conditions of obtaining 

n 

(11) (i^., i^,) = ^.-t^fc.e'(^>'=-«'=^) = 0, l<i<j<n, 

k=l 

which, after linearizing the exponential function (i.e., replacing it with its first order Taylor 
expansion) leads to (fTOl) . Therefore those phasing matrices R satisfying the real linear system 
(fTOj) lead to parametric families of complex Hadamard matrices in a neighborhood of the 
initial matrix H, up to first order; however, (11 II) is far more restrictive further decreasing 
the number of parameters in R in general. 

Note that the degree of freedom m in the defect of an n x n matrix is decreased by 2?2 — 1 
as this many parameters can always be introduced into a complex Hadamard matrix via 
multiplication by unitary diagonal matrices. This operation, however, does not yield new 
complex Hadamard matrices, up to equivalence, and therefore only dephased families are 
considered. Matrices that cannot be parametrized in any other way are called isolated. The 
most important properties of the defect are summarized in the following result from [T7t ITS]. 

Proposition 4.4. Let H be a complex Hadamard matrix. Then 

(a) d{H) is an invariant, up to equivalence, moreover d{H) = d{H*) = d{H) = d{H^); 

(b) the number of smooth parameters which can be introduced into H is at most d{H); 

(c) if d{H) = then H is isolated amongst all n x n complex Hadamard matrices. 

Note that part (c) does not require the smoothness condition and as a result it does not 
follow from part (b). 

Remark 4.5. It might be possible to use ffTOj) to extend known affine families with further 
parameters as follows. Suppose that H (a) is a dephased parametric family of complex 
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Hadamard matrices. Let > 1 be an integer and evaluate H {a) at k random points to 
obtain a series of dephased complex Hadamard matrices Hi, H2, . . . , H/^. Now consider a 
phasing matrix R (whose entries are real variables) and the resulting general parametric 
matrices Ki = Hi o EXP (iR) and the corresponding k system of linear equations given by 
ffTOj) for i = 1, . . . , k. A common solution of this k system, which should be further refined 
by ffTTl) . might result in families of complex Hadamard matrices with more free parameters 
than the initial one H {a). 

It is easy to see that a dephased complex Hadamard matrix is equivalent to a BH(g, n) 
matrix if and only if all of its entries are some gth root of unity; consequently it is easy to 
obtain all BH(g, n) matrices which are members of some parametric family H{a), stemming 
from the starting point BH(g, n) matrix H. 

5. The quaternary complex Hadamard matrices of order 10 

In this section we report on our findings regarding BH(4, 10) matrices. 

Theorem 5.1. There are exactly 10 BH(4, 10) matrices, up to equivalence, forming 7 ACT- 
equivalence classes. 

Now we present affine families containing all inequivalent BH(4,n) matrices. The first 
family, 



D'f^{a,h,c) = 
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-1 



has been obtained in [15]. It contains three out of the seven ACT-classes, namely Z)|q (1, 1, 1), 
^5^(1,1,1) and L'iJ(l,i,i). Note that its orbit can be constructed by repeatedly using 
Theorem 14.21 

The second family reads 
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1 
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whose one-parametric subfamily N[l^^{a) = N[Q^^{a, 1,1) was reported in [Ij. This matrix 

contains two ACT-classes, namely A^-[q^(1, 1, 1) and A^-[Q^(i, 1, 1), so despite the two-degree 
extension, which was obtained by a repeated application of Theorem 14.2^ no previously 
unknown BH(4, 10) matrix surfaced. 

The two remaining matrices can be obtained from complex Golay sequences [3] , and they 
belong to the family 
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The matrices G^iq{1) and G^ij{—1) are inequivalent from all previously considered exam- 
ples. We note here that the family G^iJ implies the existence of an infinite family of triplets 
of pairwise mutually unbiased bases in by a construction of Zauner [21]. 

Theorem 5.2. Each member of the ACT -equivalence classes of BH(4, 10) matrices can 
be obtained from three partially overlapping infinite affine parametric families of complex 
Hadamard matrices. These matrices, up to ACT -equivalence, can he recognized by the number 
0/3x3 vanishing minors they contain. 

It turns out that all BH(4, 10) matrices appear in the existing literature, however, the 
parametric families N^^^^ (a, b, c) and G^iq (a) are considered here for the first time. The 
authors believe that these families contain new, previously unknown BH(g, 10) matrices 
for some g > 4. Table [1] summarizes the properties of BH(4, 10) matrices for which the 
legend is as follows: the column "ACT" describes if the matrix is equivalent to its adjoint, 
conjugate or transpose, respectively, while the column "HBS" indicates if it is equivalent 
to a Hermitian matrix, if contains a n/2 x n/2 sub-Hadamard matrix (i.e., it comes from a 
doubling construction, see |12j) or if is equivalent to a symmetric matrix, respectively; the 
rest of the column headings speak for themselves. 

6. The quaternary complex Hadamard matrices of order 12 

The main result concerning BH(4, 12) matrices is the following. 

Theorem 6.1. There are exactly 319 BH(4, 12) matrices, up to equivalence, forming 167 
ACT -equivalence classes. 

Because all but four representatives of the ACT-equivalence classes contain a pair of real 
rows, almost all matrices belong to some parametric families through Lemma 14. 1[ It turned 
out that the exceptional four matrices can be parametrized as well by solving (fTOj) and 
ffTTj) . and they can be described by two parametric families (consult the families L^2a('2) 
and L^i2b{^) Appendix Therefore every BH(4, 12) matrix belongs to some parametric 
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Table 1. Summary of the BH(4, 10) matrices 



ACT Equiv. Family, ACT HBS Auto. Defect Orbit Z4 Invariant 

class classes coordinates order rank (3 x 3) 



1 


1 




YYY 


Y-Y 


64 


10 


3 


9 


2032 


2 


2,5 




NNY 


N-Y 


192 


11 


3 


9 


2496 


3 


3 




YYY 


Y-Y 


2880 


16 


3 


9 


3600 


4 


4 




YYY 


Y-Y 


32 


12 


3 


9 


2080 


5 


6,8 




NNY 


N-Y 


64 


7 


3 


9 


1568 


6 


7 




YYY 


Y-Y 


20 


8 


1 


9 


1580 


7 


9,10 




NYN 


N-N 


80 


8 


1 


9 


1600 



family. We, however, tried to minimize the number of parametric families required for the 
presentation of all BH(4, 12) matrices by attempting to describe families capturing essentially 
different properties of these matrices. We managed to describe the matrices with the aid 
of 23 families. In what follows we display four of the 23 families containing altogether a 
considerable number of distinct ACT-classes. 

The first family, stemming from the real Hadamard matrix H12, has 10 free parameters, 
and is the largest known affine family of order 12. It contains 36 distinct ACT-classes. 



Hl^^{a,b, c,d,e,f,g,h,i,j) = 
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1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


1 


1 


-1 


-1 


a 


—a 


c 


— c 


e 


e 


— e 


— e 


1 


1 


-1 


-1 


—a 


a 


c 


— c 


—e 


—e 


e 


e 


1 


1 


b 


-b 


-1 


-1 


— c 


c 


9 


-9 


h 


-h 


1 


1 


-b 


b 


-1 


-1 


— c 


c 


-g 


9_ 


-h 


h 


1 


-1 


d 


d 


-d 


-d 


cd 


—cd 


fgb 


-fgb 


-fhb 


fhb 


1 


-1 


-d 


-d 


d 


d 


—cd 


cd 


fgb 


-fgb 


-fhb 


fhb 


1 


-1 


f 


-f 


i 


—i 


—eia 


—eia 


-fgb 


fgb 


eia 


eia 


1 


-1 


f 


~f 


—i 


i 


eia 


eia 


-fgb 


fgb 


—eia 


—eia 


1 


-1 


-f 


f 


3 


-3 


eja 


eja 


-eja 


-eja 


fhb 


-fhb 


1 


-1 


-f 


f 


-j 


3 


-eja 


-eja 


eja 


eja 


fhb 


-fhb 



The second matrix has 8 parameters and contains 43 distinct ACT-equivalence classes of the 
BH(4, 12) matrices. 



Hl^'(j{a,b,c, d,e,f,g,h) = 
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1 
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—a 


— c 


-1 


-1 


-1 


-1 


1 


1 


-1 


-b 


b 


—a 


a 


-1 


e 


d 


-d 


— e 


1 


1 


-/ 


-1 


-1 


f 

J 
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J 


-/ 

J 


J 


-d 


d 


e/ 


1 


1 


f 


-1 


-1 


-f 


-f 


f 


ef 


-d 


d 


-ef 


1 


1 


— c 


b 


-6 


-1 


-1 


c 


— e 


d 


-d 


e 


1 


-1 


cgh 


bgh 


-bgh 


—agh 


agh 


-cgh 


h 


-h 


-h 


h 


1 


-1 


-cgh 


—bgh 


bgh 


agh 


—agh 


cgh 


h 


-h 


-h 


h 


1 


-1 


ch 


-bh 


bh 


-h 


-h 


—ch 


—eh 


h 


h 


eh 


1 


-1 


h 


-h 


-h 


h 


h 


h 


-h 


dh 


—dh 


-h 


1 


-1 


—ch 


h 


h 


—ah 


ah 


ch 


-h 


—dh 


dh 


-h 


1 


-1 


-h 


bh 


-bh 


ah 


—ah 


-h 


eh 


h 


h 


—eh 



The families H\ 



i'^2B ^^'^ -^12(7 iiitersect in 19 ACT-classes and they represent 60 distinct ACT- 
classes. The next family comes from Diva's general method [1] and all members of it contain 
a 6 X 6 sub-Hadamard matrix D^}\ 



, see [16]: 
D^Pia) 



D 



Dmgil,b, c,d,e,f)Di^\g) 
-Biagil, b,c,d,e J) Di^\g) 



(7) 

The family contains 20 ACT-equivalence classes, none of which is equivalent to the 
previously considered 60 classes. 

The fourth family, containing 22 ACT-classes, of which 21 is distinct from all previously 
discussed reads 



^12 {a, b,c,d,e, f,g) = 
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1 


1 
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1 
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-1 


— i 


— i 
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1 


i 


i 


— i 


— i 


ia 


ib 


-ib 


—ia 


-1 


-1 


1 


1 


-1 


-1 


e 


— e 


a 


-ib 


ib 


—a 


/ 


-/ 


1 


1 


-1 


-1 


— e 


e 


— ia 


b 


-b 


ia 




/ 


1 


1 


— i 


— i 


-1 


-1 


—a 


-b 


b 


a 


i 


i 


1 


-1 


eg 


-eg 


g 


-g 


—age 


bge 


-bge 


age 


-ifge 


ifge 


1 


-1 


eg 


-eg 


ig 


-ig 


age 


-bge 


bge 


—age 


-fge 


fge 


1 


-1 


icg 


-icg 


-ig 


ig 


dg 


-dg 


-dg 


dg 


ifge 


-ifge 


1 


-1 


-eg 


eg 


idg 


idg 


-idg 


dg 


dg 


-idg 


-dg 


-dg 


1 


-1 


-eg 


eg 


-idg 


-idg 


-dg 


idg 


idg 


-dg 


dg 


dg 


1 


-1 


-icg 


icg 


-g 


g 


idg 


-idg 


-idg 


idg 


fge 


-fge 



The four families H12B1 ^1201 -^12 ^"^^ -^12 contain 101 ACT-equivalence classes. If we 
include the adjoint, conjugate and transpose of these matrices , we end up with 186 matrices, 
up to equivalence. This is more than half of all inequivalent BH(4, 12) matrices. The further 
19 families containing the remaining BH(4, 12) matrices along with two summarizing tables 
are available in Appendix \K[ 

Theorem 6.2. Each member of the ACT-equivalence classes o/BH(4, 12) matrices can be 
obtained from 23 partially overlapping infinite affine parametric family of complex Hadamard 
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matrices. These matrices, up to ACT -equivalence, can be recognized by the number o/4 x 4 

and 5x5 vanishing minors they contain. 

7. The quaternary complex Hadamard matrices of order 14 

The exhaustive computer search yielded the following result: 

Theorem 7.1. There are precisely 752 BH(4, 14) matrices, up to equivalence, forming 298 
ACT -equivalence classes. 

Further analysis revealed 8 matrices that are isolated. For these matrices there exists a 
neighbourhood of the matrix which does not contain any inequivalent complex Hadamard 
matrices. We display here one of the 8 known isolated BH(4, 14) matrices: 



(0) 
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The defect d ( L 



(0) 



■ and hence, by Proposition I4.4[ the matrix -^[4^^ is isolated amongst 
all 14 X 14 complex Hadamard matrices. Three additional inequivalent isolated matrices can 
be obtained by considering the adjoint, conjugate and transpose of L 



(0) 
14A- 



Theorem 7.2. There are at least 8 isolated BH(4, 14) matrices, up to equivalence, forming 
two ACT -equivalence classes. 

As we have found isolated matrices it follows that it is not possible to come up with a 
universal parametrization scheme for the BH(4, n) matrices. This stands in contrast to the 
real Hadamard matrices which can be parametrized always (cf. Lemma [4. ip . 

Theorem 7.3. All BH(4, 14) matrices, up to ACT -equivalence, can be recognized by the 
number o/4 x 4 vanishing minors they contain. 



Although by using Theorem 14. 2 1 we were able to obtain various parametric families starting 
from many of the BH(4, 14) matrices we constructed, we could not introduce affine parame- 
ters into matrices having relatively small defect. Whether or not they are isolated remains 
an open problem. 

For an illustration of the parametrization of BH(4, 14) matrices, we display here a six- 
parameter family of complex Hadamard matrices: 



(a, b, c, d, e. 



/) 
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The starting-point matrix D[J{1, 1, 1, 1, 1, 1) is a symmetric conference matrix, which was 
considered along with the five-parameter subfamily Df^{a,b,c,d,e) = Df^{a,b,c,d,e,l) 
in [15] . The extra parameter was found by the argument outlined in Remark 14.51 we have 
evaluated the known five-parameter family D^^ (a, b, c, d, e, 1) at random fourth roots of unity 
to obtain some BH(4, 14) matrices, say Hi, H2, . . . , Hiq, and considered the corresponding 10 
instances of f lTU]) featuring the same phasing matrix R in all ten cases. These extra equations 
heavily constrained the matrix R resulting in the one-parametric extension we discovered. 
Despite its large degree of freedom the family D^'^ contains 14 ACT-classes only, and it seems 
that a given affine parametric family cannot contain significantly more different ACT-classes 
in order 14. 

For a summarizing table highlighting the main features of BH(4, 14) matrices and for 
additional examples of parametric families consult a Web site [9j. All results of this work 
are available in full detail at the Web site which describes all inequivalent BH(4, n) matrices, 
parametric families and summarizing tables for 10 < n < 14. 
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Appendix A. A list of parametric families of BH(4, 12) matrices 

In the following we display the 19 parametric families of order 12 which were not shown 
in Section O These 19 families with the families h[2b, h[^^,, ^ (defined in Section 
in]) account for all BH(4, 12) matrices, up to ACT-equivalence: 
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Table H] summarizes the relations between BH(4, 12) matrices and parametric families by 
showing which ACT-classes are included in each of the 23 parametric families. An ACT- 
class number is underlined if the ACT-class belongs to only one of the parametric families 
listed in the table. 

Table |3] summarizes the properties of BH(4, 12) matrices. 
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Family ACT-classcs 

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 21, 

22, 23, 28, 29, 45, 46, 47, 48, 49, 59, 62, 65, 66, 70, 78, 79, 81 
h[1^^ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 

27, 28, 32, 33, 34, 37, 38, 39, 40, 50, 51, 52, 53, 55, 57, 58, 60, 61, 63, 64 
D^j' 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 

113, 114, 115, 117, 127, 128, 129, 130, 131, 138 
X[l^ 44, 75, 76, 107, 116, 118, 120, 121, 122, 123, 124, 

125, 126, 133, 134, 135, 136, 137, 139, 140, 154, 160 
l[2a 164. 165 
-^1213 166.167 

5*2 4 30,31,68,69,74,87,88,156,157 

^iJb 103, 104, 105, 106, 107, 108, 116, 117, 121, 122, 123, 124, 125, 126, 

127, 129, 130, 132, 133, 134, 141, 142, 143, 144. 153. 158. 161, 162. 163 

sjjc 32, 33, 35, 36, 37, 38, 58, 60, 61, 63, 64, 71, 72, 80 

Sjjn 1,2,3,6,10,21,24,41,56,73 

1, 2, 3, 4, 6, 10, 15, 21, 24, 26, 39, 40, 41, 42, 51, 67 

B^2F 1, 2, 3, 8, 21, 22, 23, 24, 25, 37, 38, 74, 83, 93, 97, 111, 138 

B^2G 1, 2, 3, 8, 9, 21, 22, 23, 24, 25, 34, 37, 38, 73, 84, 85, 94, 139, 140 

B[f^ 1,2,3,8,9,34,86,147 

Bjjr 30.31.89.91,148, 150 

b\Yj 13, 35, 47, 59, 78, 81, 92, 95, 105, 121 

^12K 13, 36, 47, 59, 78, 81, 90, 96, 108, 123 

b[^^^ 1,2,3, 4, 21, 22, 23, 26, 32, 33, 58, 98, 99, 100, 101, 141, 142, 143, 144 
■^1 2M 30, 43, 87, 89, 102, 118, 119, 133, 145 
sjjL 1 19, 145, 146, 147, 154 
BiJo 148, 149 

B[fp 74, 75, 93, 107, 109, 111, 112, 122, 151, 152, 156, 157 
^120 30. 43. 54, 68. 77, 82, 118, 119. 155. 159 
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NYN 


NNN 


32 


11 


7 


5 


(27878 


34616) 


95 


149, 294 


b'*' (1,1,-1,1) 

12 J V ' ' ' / 


NYN 


NNN 


32 


11 


4 


6 


(20998 


31464) 


96 


150, 291 


B(tt(l, 1,1,1) 

12k V ' ' ' / 


NYN 


NYN 


32 


13 


4 


6 


(21566 


46832) 


97 


151, 273 


-Bioi-'fi, —1, i, — i, —1, i, — i) 


NYN 


NYN 


160 


17 


7 


5 


(30350 


60480) 


98 


152, 293 


si~l(-l,-l,-l,i, -i) 

12Z/ ' ' ' ' / 


NYN 


NYN 


32 


11 


5 


6 


(21118 


36288) 


99 


153, 275 


BrJr (1, 1, -1, 1, -1) 

12L V ' ' ' ' ■' 


NYN 


NYN 


32 


9 


5 


6 


(21478 


29888) 


100 


154, 276 


sJSl (1, 1, 1, 1, 1) 


NYN 


NNN 


32 


9 


5 


6 


(21070 


28800) 


101 


155, 292 


Sj^^r (1,1,-1, -i,-i) 

12Z/ V ' ' ' ' / 


NYN 


NNN 


32 


13 


5 


6 


(21118 


36864) 


102 


156,312 


b\1\.(i, 1, 1) 


NYN 


NYN 


4 


11 


3 


7 


(18750 


31342) 


103 


157 


d\V(i, -1, -1, 1, -1, -1, 1) 

]^2 \ ' ' ' ' ' ' / 


YYY 


YYY 


960 


27 


7 


6 


(33270, 


106080) 


104 


158, 280 


D(^>(-l,-l,l,i, -1,-1,-1) 


NYN 


NYN 


160 


17 


7 


5 


(29510 


85680) 


105 


159,225 


D$^'(l,-i,-l,i, -1,1,-1) 


NYN 


NYN 


64 


19 


7 


6 


(27550 


69280) 


106 


160, 284 


£)JJ(l,-l,-l,l,-l,-i,-l) 


NYN 


NYN 


32 


19 


7 


5 


(29078 


70144) 


107 


161, 183 


d(^)(1, 1,1, 1,-1, 1,-1) 


NYN 


NYN 


192 


25 


7 


6 


(32622 


79488) 


108 


162, 252 


-^1? (-1, -1, -i, 1, -i, -1, -1) 


NYN 


NYN 


64 


19 


7 


6 


(27558 


74208) 


109 


163,216 


D<J(i,i,l,l,l,-i,-l) 


NYN 


NYN 


16 


19 


7 


6 


(27442 


66528) 


110 


164, 286 


D(^)(i,i,l,i,i,l,-1) 


NYN 


NYN 


32 


17 


7 


5 


(28982 


67952) 
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ACT 




Pa milv 

J. ClilllllJ'j 


ACT 


HBS 


Auto. 


Dpfprt 


Orbit 


Z4 


Invariant 


class 


classes 


coordinates 






order 






rank 


(4x4 


5x5) 


111 


165, 218 


D'i'^{-1, -i, 1, -1, 1, i,i) 


NYN 


NYN 


192 


25 


7 


6 


(32358 


73632) 


112 


166,219 


D\J (i, i, i, 1, i, i, — i) 


NYN 


NYN 


16 


17 


7 


6 


(27406, 


66632) 


113 


167, 303 


Di^^(l,l, -l,-i,i,i, -1) 


NYN 


NYN 


32 


13 


7 


5 


(27070, 


67008) 


114 


168, 304 


D5?(i, -1,1,1,1,1, -1) 


NYN 


NYN 


48 


13 


7 


5 


(26946, 


64392) 


115 


169, 307 


d\V (-1,1, -1, -1, 1, 1, -1) 


NYN 


NYN 


96 


13 


7 


5 


(26958, 


70512) 


116 


170, 196 


X$?(-l, -1,1, -1,1,1, 1) 


NYN 


NYN 


64 


19 


7 


6 


(26870, 


50752) 


117 


171, 186 


(7") 

D\n (i- i- — i, —1, — 1) 

\2 V 7 7 ) 7 7 7 / 


NYN 


NYN 


64 


19 


7 


6 


(27646, 


76160) 


118 


172, 197 


xf? (1,-1, 1,1, -1,1,-1) 


NYN 


NNN 


16 


15 


7 


5 


(27314, 


51064) 


119 


173, 212 


12M V 5 1 y 


NYN 


NNN 


16 


13 


3 


5 


(27056, 


43288) 


120 


175, 208, 267, 274 


xiP(l, -1, -1, 1, 1, 1, 1) 


NNN 


NNN 


32 


9 


7 


5 


(25826, 


42248) 


121 


176, 253 


x[l\i, 1, -1, 1, 1, -1, 1) 


NYN 


NYN 


64 


17 


7 


6 


(27142, 


47472) 


122 


178, 269 


Xj2'(l,l, -1.1, -1.-1,-1) 


NYN 


NYN 


32 


19 


7 


5 


(28254, 


58488) 


123 


179, 224 


X)n (i, — i, 1, — 1, i) 


NYN 


NYN 


64 


17 


7 


6 


(26350, 


45024) 


124 


181, 266 


X^V (1,1,-1,1,-1,-1,1) 


NYN 


NYN 


32 


17 


7 


5 


(27678, 


45536) 


125 


182 


X(^'(l, -1,-1, -1,-1, -1,-1) 


YYY 


YYY 


192 


27 


7 


6 


(32406, 


57120) 


126 


184, 199 


xi? (1, 1, -1,-1,-1, -1, 1) 


NYN 


NYN 


64 


19 


7 


6 


(26814, 


47424) 


127 


185, 190 


D^Vii, 1, -1, 1, -1, -1, -1) 


NYN 


NYN 


64 


19 


7 


6 


(27494, 


71648) 


128 


187, 282 


D^Vih -1, 1, -1, 1, -1, 1) 


NYN 


NYN 


160 


17 


7 


5 


(29030, 


74720) 


129 


188, 300 


Di? (-1,1, 1,-1, -1,-1,1) 

]^2v ' 7 5 7 5 5 / 


NYN 


NYN 


32 


17 


7 


5 


(27078, 


70272) 


130 


189, 306 


DjI'd, -1,-1,1, 1,1,1) 


NYN 


NYN 


96 


21 


7 


5 


(27246, 


77088) 


131 


191,302 


D$?(l, -1,-1, 1,1,-1.-1) 


NYN 


NYN 


48 


21 


7 


5 


(26814 


66000) 


132 


192, 193 


b$od(i, -1, -1, -1, 1) 

12H \ ' ' ' ' / 

x,v(i,i, -1,-1,1, 1,-1) 

\ } 7 } 5 3 7 / 


YNN 


YYN 


16 


13 


5 


8 


(20454, 


40556) 


133 


195, 305 


NYN 


NYN 


32 


17 


7 


5 


(26090, 


51248) 


134 


198, 301 


(1,1, -1,-1, -1,1,1) 

]_2 \ ' ' ■ ' 7 5 / 


NYN 


NYN 


32 


19 


7 


5 


(26006 


46256) 


135 


200, 217 


xi? (-1, 1, -1, -1, -1, 1, -1) 

12 V -^)-*5 -^5 -^5 -^J-^) 


NYN 


NNN 


16 


13 


7 


5 


(27126, 


44376) 


136 


201, 285 


xj?(l, -1, 1,1, 1, -1,1) 


NYN 


NNN 


32 


11 


7 


5 


(27390, 


40048) 


137 


202, 281 


x(?(-l, 1, 1, -1, -1, -1, -1) 

\2, V 7 ' 5 7 ' 7 / 


NYN 


NNN 


32 


11 


7 


5 


(27614, 


45968) 


138 


203, 258 


d(?(1, -1,1,-1, -1,1, 1) 

]_2 \ 7 11 1 1 1 J 


NYN 


NYN 


960 


25 


7 


6 


(31950, 


76800) 


139 


204, 257 


x'j'ci, -1,1,-1,-1,-1,1) 


NYN 


NNN 


192 


19 


7 


6 


(31302, 


45600) 


140 


205, 242 


X<|' (-1,-1, -1,1, -1,-1,1) 


NYN 


NNN 


192 


19 


7 


6 


(30990 


41088) 


141 


206, 262 


12s 


NYN 


NYN 


64 


15 


5 


6 


(20318 


42192) 


142 


207, 259 


5(5) (_j^ _^ _^ _^ 


NYN 


NYN 


64 


17 


5 


6 


(20062 


38544) 


143 


210,261 


sioD(l, 1,1,1, 1) 

12o V ' ' ' ' / 


NYN 


NYN 


64 


15 


5 


6 


(20606, 


37936) 


144 


211,260 


Slond, -1, -1, -1, 1) 

12s \ ' ' ' ' -' 


NYN 


NYN 


64 


17 


5 


6 


(20654, 


36592) 


145 


213, 308 


B!oL(i, 1, -1) 

12M \ ' ' / 


NYN 


NYN 


48 


19 


3 


5 


(25542 


42240) 


146 


214, 283 


5$^) (_1 1 i i 1 1) 
12iV V ' ' ' ' ' / 


NYN 


NNN 


40 


9 


6 


5 


(27260, 


37940) 


147 


215, 256 


s!oL(-i, -1,1,1,-1,1) 

12h V ' ' ' ' ' / 

BllUi, 1,1,1) 

IzJ \ ' ' ' / 


NYN 


NNN 


240 


17 


6 


5 


(31140, 


37080) 


148 


220, 223 


NYN 


NYN 


16 


15 


4 


8 


(19382 


29480) 


149 


221,311 


B$2o(l,l>l) 


NYN 


NNN 


16 


11 


3 


7 


(18036, 


22832) 


150 


222 


1.1.1) 


YYY 


YYY 


16 


19 


4 


8 


(19982, 


33920) 


151 


226,227,287,289 


B^^^p(l, 1,1,1) 


NNN 


NYN 


16 


10 


4 


6 


(21018 


41104) 


152 


235, 290 


Bj^' (1,1,1,-1) 


NYN 


NYN 


8 


15 


4 


6 


(21182, 


40566) 


153 


236, 313 


sjf (-1,-1,1,1,1) 


NYN 


NYN 


8 


11 


5 


7 


(18600, 


36956) 


154 


237, 248 


x|2^(-l, 1,-1, 1,1, -1,1) 


NYN 


NNN 


48 


15 


7 


5 


(29142 


49512) 


155 


238, 246 


B^2q(1. 1.1.1) 


NYN 


NNN 


16 


11 


4 


5 


(24692, 


38312) 


156 


239, 245 


Si2a(!. 1.-5.1) 


NYN 


NNN 


16 


11 


4 


7 


(19578, 


31544) 


157 


240, 244 


5^*^(1,-1,-1,1) 
(1,-1, -1,1,1) 


NYN 


NNN 


16 


11 


4 


7 


(19794 


30180) 


158 


241, 298 


NYN 


NYN 


32 


11 


5 


6 


(20686, 


34552) 


159 


247, 250 


b|2q(i.i.-i.i) 


NNY 


NNY 


32 


11 


4 


5 


(24960, 


45408) 


160 


249 


xJj'C-l.i. 1. -l.l.-l.i) 


YYY 


NNY 


48 


15 


7 


5 


(29364 


56832) 


161 


251, 296 


^I2s(-l.-l.l.-l.i) 


NYN 


NYN 


32 


13 


5 


6 


(21222, 


49720) 


162 


254, 297 


(-1,1,1,1,1) 


NYN 


NYN 


32 


13 


5 


6 


(21806, 


46744) 


163 


255, 299 


5(5^^(1,1,-1,1,-1) 


NYN 


NYN 


32 


11 


5 


6 


(21334 


32856) 


164 


314,317 


NYN 


NYN 


8 


7 


1 


9 


(14130 


19680) 


165 


315, 319 




NYN 


NNN 


24 


9 


1 


9 


(14016, 


17172) 


166 


316 




YYY 


NYY 


8 


11 


1 


9 


(14284 


23940) 


167 


318 




YYY 


NNY 


24 


9 


1 


9 


(13440, 8640) 



